Classical realization of two-site Fermi-Hubbard systems 
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Quantum-classical analogies have been explored on 
many occasions to mimic and visualize in a purely clas- 
sical setting the dynamical aspects embodied in a wide 
variety of quantum systems ^ [2]. In particular, in the 
past two decades engineered photonic lattices have pro- 
vided a useful model system to investigate wave optics 
analogous of solid state phenomena [2H5]. Most of the 
optical analogues of solid-state phenomena observed so 
far, including electronic Bloch oscillations [3| [6], Zener 
tunneling [7:, dynamic localization [5], Anderson local- 
ization [9 1, Rabi flopping [TO], and topological photonic 
crystals refer to single-particle phenomena and are 
based on the formal similarity between the paraxial op- 
tical wave equation in photonic lattices and the nonrela- 
tivistic Schrodinger equation of a single particle in peri- 
odic potentials 2J. However, much of the richer physics 
in condensed-matter comes from many-body phenomena 
and electron correlations. The simplest and paradigmatic 
model which describes correlation effects of electrons in a 
lattice, arising from the the competition among chemical 
bonding, Coulomb repulsion and Pauli exclusion prin- 
ciple, is perhaps provided by the Hubbard model (HM) 
|12j . This model is capable of capturing some many-body 
aspects of the electronic properties of condensed matter, 
such as metal-insulator transitions, itinerant magnetism, 
and electronic superconductivity (see, e.g., [12] [T3] and 
reference therein) . In spite of the simplicity of its Hamil- 
tonian structure, very few exact results are known for 
the HM, mainly for finite clusters or for the infinite one- 
dimensional chain [T3T - fT5] . The simplest solvable and 
nontrivial system, which can still capture some of the 
main relevant properties of larger clusters and of the in- 
finite chain, is provided by the two-site Hubbard Hamil- 
tonian (see, for instance, [H]). The two-site HM, being 
exactly solvable, has been considered by several authors 
as a simplified theoretical model [Tt>H2"2"] . In particular, 
it is useful as a toy model for understanding the bind- 
ing of molecules like H2 [T7Tll9j . and it was proposed to 
model electron-molecular vibration coupling in organic 
charge-transfer salts J2U] and the electronic structure in 
7r systems [3T]. Since photons are bosons and they do 
not interact when propagating in linear optical struc- 
tures, one would expect that photonics is not a suited 
system to simulate in a classical setting the physics of in- 



teracting electrons in solids. In recent works [23], it has 
been pointed out that photonic structures could provide 
a noteworthy laboratory system to simulate the physics 
of few interacting bosons in the framework of the Bose- 
Hubbard model. In this Brief Report it is shown that 
light transport in suitably engineered coupled waveg- 
uide structures can mimic the dynamics of interacting 
fermions as well. In particular, an optical realization of 
the two-site HM is proposed, in which light propagation 
in four evanescently-coupled waveguides reproduces the 
temporal dynamics of the occupation number amplitudes 
of the electrons in the two-site potential. 

The HM describing electron dynamics in a one- 
dimensional chain of N potential sites with nearest neigh- 
boring hopping is defined by the Hamiltonian (see, for 
instance, [14]) 

JV-l N 
3=1 <r=t4 3=1 

(1) 

where k is the hopping amplitude between adjacent sites, 
U is the on-site Coulomb interaction strength, a}- a is 
the fermionic creation operator that creates one elec- 
tron at site j with spin a (j = 1, 2, ...,N, a =^.\), and 
"hj,cr = o,j o-%,ct are the particle number operators. The 
fermionic operators a}- a and fi Ji<T satisfy the usual anti- 
commutation relations {fit ff , fi|. } = {&j tCr , ak, P } = and 

{dj^,a\ p } = Sj t kO~a, P - The space of states of the HM is 
spanned by all linear combinations of Wannier states of 
the form pH] 

\m,n 2 , ...,njv,mi,m2,...,mjv) = |n, m) = (2) 

— "l,l U 24 ■■■ u N,i "l,t "2,t •■■ U JV,T I ' 

where |0) is the vacuum state. The state |n, m) corre- 
sponds to rij electrons occupying the site j with spin ! 
and m,j electrons occupying the site j with spin f . Owing 
to the anti-commutation rules of the Fermi operators, the 
integers n 3 - and m 3 - can take only the two values and 
1, according to the Pauli exclusion principle. Hence, the 
number of all linearly independent Wannier states is 2 2N . 
If the state vector \ijj(t)) of the system is decomposed on 
the Wannier basis, \tp(t)) — m /(n, m, t)|n, m), the 
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FIG. 1. (a) Wannier basis of the two-site Hubbard Hamil- 
tonian corresponding to two electrons with opposite spins, 
and (b) corresponding optical realization based on four 
evanescently-coupled optical waveguides. The coupling rate 
between nearest waveguides is k, whereas the propagation 
constant of the modes of waveguides 1 and 2 is shifted by U 
from the one of waveguides 3 and 4. The waveguides 1,2,3 and 
4 correspond to the four quantum states 1 1, 0, 1,0), |0, 1,0, 1), 
|0, 1,1,0) and |1, 0,0,1), respectively. 



evolution equations for the 2 2N occupation amplitudes 
/(n, m, i) are formally given by (assuming H = 1) 



. G?/(n, m, t) 
dt 



s,q 



n,m|#|s,q)/(s,q,i). 



(3) 



Since the total number of electrons N t and total num- 
ber of electrons with spin f C^f) an d J, (-^Yj.) are con- 
served quantities for the Hubbard Hamiltonian [53], the 
amplitude equations (3) are decoupled into a set of equa- 
tions acting on different sub-spaces of the Hilbert space. 
Each sub-space is defined by the Wannier states with 
an assigned number of electron Nt = iVj, + Nf, with 
N-f electrons with spin f and N± electrons with spin J,; 
the number of amplitudes in such a subspace is hence 

' N \ ( N \ 

\ The two-site HM [i.e. TV = 2 in Eq.(l)] 



provides the simplest and exactly solvable model which 
can still capture some of the main relevant properties of 
larger clusters and of the infinite chain. The two-site 
HM has been investigated by several authors [T6H22] and 
proposed as a toy model for understanding the binding 
of molecules like H2 [TTHHj] , to model electron- molecular 
vibration coupling in organic charge-transfer salts [2TJ] . 
and to describe the electronic structure in ir systems 
[2"T] . Here we propose an optical realization of the two- 
site Hubbard Hamiltonian based on light transport in 
evanescently-coupled optical waveguides which is capable 
of mimicking the temporal dynamics of the Hubbard sys- 
tem in the Wannier basis representation (3). For the two- 
site Hubbard Hamiltonian, there arc 16 possible config- 
urations of electrons on two sites: one with no electrons, 
four with one electron (an up or down electron on each 
of the two sites) , six with two electrons (one with two up 
electrons on different sites, one with two down electrons, 
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FIG. 2. Temporal evolution of (a) the return probability 
P(t), and (b) the spin imbalance Ni2(t) for the two-site Hub- 
bard model for U/k = 0.5. The system is initially prepared 
in the state |0, 1, 1, 0). 
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FIG. 3. Same as Fig.2, but for U/k, = 5. 



and four with an up electron and a down electron), four 
with three electrons, and one with four electrons. The 
most interesting dynamics is provided by the sub-space 
consisting of two electrons with opposite spins, i.e. to 
N t = 2 and N^ — = 1, which is spanned by the four 
states |1, 0,1,0), |0, 1,0,1), |0, 1,1,0) and 1 1,0, 0,1) with 
amplitudes c x {t) = f(l, 0, 1, 0, t), c 2 (t) = /(0, 1, 0, 1, t), 
c 3 (t) = /(0, 1,1,0,*) and c 4 (t) = /(l,0,0,l,i) [see 
Fig. 1(a)]. In this case, the coupled equations (3) for the 
amplitudes q(£) (/ = 1,2,3,4) read explicitly 
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An optical realization of the Hamiltonian system 
(4) is provided by propagation of monochromatic light 
waves at wavelength A in four evanescently-coupled op- 
tical waveguides in the geometrical setting shown in 
Fig. 1(b). In the optical structure, the time t repre- 
sents the spatial propagation distance along the waveg- 
uide axis, whereas the Wannier amplitude q corre- 
sponds to the modal amplitude of light trapped in the 
l-th waveguide (I = 1,2,3,4). In fact, in the tight- 
binding approximation the spatial part of the electric 
field £{x, y, t) of the optical wave propagating along the t 
spatial direction of the guiding structure can be written 
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as £(x,y,t) ~ Xw=i c i(t) u i( x : V) exp(2ninot / X) , where 
(x, y) are the spatial coordinates transverse to the op- 
tical t axis, ui(x,y) is the spatial modal profile of the 
Z-th waveguide, and hq is an effective reference mode in- 
dex. The spatial evolution of the modal amplitudes q 
arising from the weak evanescent coupling of adjacent 
waveguides is governed by coupled mode equations 
which are formally analogous to Eqs.(4), provided that 
the cross-coupling between waveguides 1 and 2, and be- 
tween waveguides 3 and 4 in Fig. 1(b), is negligible. In 
the optical setting, the hopping rate k entering in Eq. (4) 
is analogous to the spatial tunneling rate of light waves 
between two adjacent waveguides arising from evanes- 
cent field coupling, whereas the on-site Coulomb inter- 
action strength U corresponds to a shift of the propa- 
gation constants of waveguides 1 and 2 as compared to 
waveguides 3 and 4 [see Fig. 1(b)]. The optical structure 
shown of Fig. 1(b) could be easily realized in fused silica 
by the recently-developed femtosecond writing technique, 
in which the propagation constant shift U is realized by 
varying the writing speed of waveguides (see, for instance, 
0). 

The energies of the two-site Hubbard Hamiltonian 
can be calculated analytically as the eigenvalues of the 
4x4 matrix entering in Eq.(4), and read explicitly (see, 
for instance, [20]) E 1 = 0, E 2 = U, E 3 = (U/2) + 
y/(U/2) 2 +4k 2 , and E 4 = (U/2) - y/(U/2) 2 + 4k 2 . In 
the optical analogy, such energies correspond to the prop- 
agation constant mismatch of the various supermodes of 
the coupled waveguides. It is worth noticing that the en- 
ergy spectrum of the simple two-site HM contains some 
important physical features of more complex chains, such 
as the onset of metal-insulator transition in a half-filled 
linear chain as the parameter k/U is varied [16| . To see 
this, let us consider the limit of small k/U, and let us 
expand the sector of Hilbert space to include all sec- 
tors with two electrons by adding the states 1 0011) and 
1 1100). These states are eigenstates of H with eigen- 
value 0. All together, the two electron space of the 
two site HM has four 'small' eigenvalues 0, 0, and 
(U/2) - ^(U/2) 2 + Ak 2 ~ -Ak 2 /U, and two 'large' ones 
U and (U/2) + ^/(U/2) 2 + Ak 2 ~ U. The large eigen- 
values are associated with eigenvectors whose compo- 
nents have significant mixtures of the states with dou- 
bly occupied sites. The existence of the two groups of 
states whose eigenvalues are separated by U is a reflec- 
tion of the upper and lower Hubbard bands in a lat- 
tice. The 'Mott-Hubbard' gap in the spectrum gives rise 
to a metal- insulator transition [14] . For a given initial 
state \ip(t — 0)) = \ipo), there are two interesting observ- 
ables related to the dynamical evolution of the two-site 
HM, namely the return probability P(t) — \(ip(t)\ipo)\ 2 
and the spin imbalance between the two sites, Ni 2 (t) = 
(l/2)(ip(t)\ni : f — ni : i + n 2 ,i — h2,f\ip(t)) (see, for instance, 
[22]). The latter describes the exchange dynamics of the 
two spins f and I, located at the two sites. As an ex- 
ample, Figs. 2 and 3 show a typical behavior of the re- 
turn probability P(t) and spin imbalance N 12(f) for the 



two-site HM corresponding to a weak (U/k — 0.5, Fig. 2) 
and a strong (U/k = 5, Fig. 3) on-site Coulomb interac- 
tion for a system initially prepared in the Wannier state 
|0, 1,1,0), i.e. for the initial condition c;(0) = <5;.3. In 
the Wannier basis representation, the return probability 
P(t) and spin imbalance Ni 2 (t) take the simple form 

4 

P(t) =$>?(0)Q(i)| 2 , N 12 (t) = | C3 (<)| 2 -|c 4 W| 2 . (5) 
1=1 

Note that in the optical analogue the spin imbalance 
has a very simple meaning: it is just the power imbal- 
ance of light between waveguides 3 and 4. Moreover, if 
the system is initially prepared in one of the Wannier 
state, the return probability is simply mapped into the 
fractional optical power trapped in the initially excited 
waveguide. For example, if the system is initially pre- 
pared in the state |0, 1, 1, 0) as in Figs. 2 and 3, one has 
P(t) = I C3 (t) 1 2 . It is worth noticing that, the optical ana- 
logue of the strong on-site Coulomb interaction regime 
U/k ^> 1 (as in Fig. 3) corresponds to a nearly-sinusoidal 
exchange of optical power between waveguides 3 and 4, 
like in an ordinary synchronous optical direction coupler 
|25) . In fact, in the large U/k limit and for the initial 
condition cx(0) = c 2 (0) = 0, the amplitudes c\ and c 2 re- 
main small and can be eliminated from the dynamics by 
standard perturbation methods. This yields the reduced 
dynamical equations for the amplitudes c 3 and c 4 

i ~n~ — ~«eC4 + Sc 3 (6) 

at 

1— ~ -K e c 3 + dc 4 (7) 
at 

where n e = 2k 2 /U is an effective coupling constant 
and (5 = — K e a common propagation constant detun- 
ing. Hence light coupling in the four-waveguide structure 
in the large on-site Coulomb interaction regime is like 
the one of an ordinary two-waveguide directional cou- 
pler, showing Rabi-like exchange of the optical power 
between the (not directly coupled) waveguides 3 and 4 
mediated by the two off-resonance waveguides 1 and 2. 
The observed oscillatory power exchange between two 
uncoupled waveguides, arising from indirect coupling via 
weakly excited off-resonance waveguides, is analogous to 
the process of two-photon Rabi oscillations observed in 
Ref.[H]. 

In conclusion, a simple classical realization of the two- 
site Fermi- Hubbard Hamiltonian, based on light propaga- 
tion in evanescently-coupled optical waveguides, has been 
theoretically proposed. While previous theoretical and 
experimental studies on optical simulations of quantum 
phenomena in solid-state physics have been concerned 
with single-particle phenomena, here it has been shown 
that photonics can provide a laboratory tool to visual- 
ize and simulate in simple optical settings the dynamical 
aspects embodied in the physics of interacting fermionic 
systems. The present study has been focused on a sim- 
ple Hubbard system, however it is envisaged that pho- 
tonic simulators of others and more complex models of 
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interacting fermions can be realized. Possible extensions 
of the present study include the photonic realization of 
the Hubbard-Holstcin Hamiltonian [37], which describes 
bipolarons dynamics arising from electron-phonon cou- 
pling, and the Hubbard- Anderson Hamiltonian [35] de- 
scribing the dynamics of two electrons on a linear chain 
with long-range correlated disorder and on-site Hubbard 
interaction. The Hubbard- Anderson Hamiltonian can be 
simulated using a two-dimensional square array of waveg- 
uides with engineered propagation constants, which can 
be used to test in an optical setting the interplay be- 



tween disorder, localization and electron-electron interac- 
tion. Phonon-electron coupling dynamics for a two-site 
Hubbard-Holstein Hamiltonian can be realized by cou- 
pling two of the four waveguides of Fig. 1(b) with two 
semi-infinite linear arrays with non-homogeneous hop- 
ping rates, which simulate the vibrational (phonon) de- 
grees of freedoms of the two sites. 
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